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In a recent paper [I. P. Neupane and D. L. Wiltshire, Phys. Lett. B 619, 201 (2005).] we have 
found a new class of accelerating cosmologies arising from a time-dependent compactification of 
classical supergravity on product spaces that include one or more geometric twists along with non- 
trivial curved internal spaces. With such effects, a scalar potential can have a local minimum with 
positive vacuum energy. The existence of such a minimum generically predicts a period of accelerated 
expansion in the four-dimensional Einstein-conformal frame. Here we extend our knowledge of these 
cosmological solutions by presenting new examples and discuss the properties of the solutions in a 
more general setting. We also relate the known (asymptotic) solutions for multi-scalar fields with 
exponential potentials to the accelerating solutions arising from simple (or twisted) product spaces 
for internal manifolds. 
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I. INTRODUCTION 



o 

The possibility that fundamental scalar helds which are uniform in space play a preeminent role on cosmological 
scales has been confirmed by a decade of observations. Most recently the WMAP measurements of fine details 
J> . of the power spectrum of cosmic microwave background anisotropies 0] have lent strong support to the idea that 
the universe underwent an early inflationary expansion at high energy scales. The WMAP data, along with the 
independent observations of the dimming of type la supernovae in distant galaxies 0] are also usually interpreted as 
an indication that the universe is undergoing accelerated expansion at the present epoch, albeit at a vastly lower rate. 

As in the case of the much earlier period of early universe inflation, the most natural explanation for the repulsive 
force responsible for accelerating cosmologies would be a fundamental vacuum energy, possibly in the form of one 
or more dynamical homogeneous isotropic scalar fields. The nature of this dark energy, which constitutes of order 
j 70% of the matter-energy content of the Universe at the present epoch, constitutes a mystery whose explanation is 
-<-j J . possibly the greatest challenge faced by the current generation of cosmologists. 

Qh! In the ACDM model, the dark energy at the present epoch is attributed purely to a constant vacuum energy (or 
cosmological constant), or equivalently a homogeneous isotropic fluid whose pressure, p, and energy density, p, are 
related by p = —p = A. This is only the simplest (and perhaps most common) explanation. From a field theoretic 
viewpoint it would be perhaps more natural to attribute the dark energy to one or more dynamical scalar fields 0, Q . 
Many such quintessence models have been studied, with scalar potentials which range from completely ad hoc ones to 
those with various theoretical motivations. Such motivations are often more than not phenomenological. For example, 
ultra light pseudo Nambu-Goldstone bosons can give realistic cosmologies [jg, even if one does not specify exactly 
what fundamental theory such scalars belong to. 

Fundamental scalar fields are of course abundant in higher-dimensional theories of gravity. The typical scalar 
potentials that one obtains by dimensional reduction, namely exponential potentials have been widely studied Q but 
generally without regard to the restrictions on the sign of the potential and magnitude of the coupling constants that 
arise from compactifications or the higher-dimensional geometry. In the case of multiple scalar fields, for example, 
attention has been focused on simplified potentials in which each scalar appears only once 0, 13, EJ El El • 

A number of models with a fundamental higher-dimensional origin, which accommodate a 4-dimensional universe 
with accelerating expansion have been studied over the past two decades 0, El El , usually in relation to inflationary 
epochs in the very early universe. A general feature of these models was the presence of additional degrees of freedom, 
such as fluxes or a cosmological term. Until relatively recently it was not believed that one could obtain accelerating 
universes in the Einstein frame in four dimensions by compactification of pure Einstein gravity in higher dimensions 
[l5j . In fact, the result was elevated to the status of a no-go theorem 
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Townsend and Wohlfarth 17] circumvented the no-go theorem by relaxing its conditions to obtain time-dependent 
(cosmological) compactifications ol pure Einstein gravity in arbitrary dimensions, with negatively curved internal 
space, which are made compact by topological identifications. Their exact solutions exhibit a transient epoch of 
acceleration between decelerating epochs at early and late times, albeit with the number of e-folds of order unity. 
Applied to late-time cosmologies, this could mean that our present cosmic acceleration is not eternal but switches off 
in a natural way. 

Transient acceleration was subsequently shown to be a generic feature of many supergravity compactifications 

El E El El El El in ei ei ei ei , some models including higher dimensional fluxes appropriate to particular 

supergravity models, and others with more complicated internal product spaces with negatively curved factors. The 
generic nature of transient acceleration is understood also for some simple potentials arising from group manifold 
reductions m , and a compactification of M-theory on a singular Calabi-Yau space |3(| ■ All these solutions typically 
exhibit a short period of accelerated expansion. Other approaches to accelerating cosmologies involve j3l| a compact- 
ification of string/M theory with higher order curvature corrections, such as (Weyl) 4 terms in type II string theory, 
see also a review |32l . In this rather complicated scenario, the analysis is often restricted to asymptotic solutions of 
the evolution equations, which involve a fine-tuning of the string coupling or the Planck scale in the higher dimensions 
(10 or 11). It is not immediately clear if any asymptotic solutions would be available within the full string theories. 

In a recent paper (3j| we have found a new class of supergravity cosmologies arising from time-dependent compact- 
ifications of pure Einstein gravity, whose existence is somewhat counter to the intuition of Ref. [l7|. In particular, we 
have exhibited exact solutions which circumvent the no-go theorem |l6j . while retaining Ricci-flat internal spaces. 
This is possible through the introduction of one or more geometric twists in the compact space. The observation 
that acceleration is possible even for Ricci-flat cosm olog ical compactifications is not new, which only required the 
introduction of external fluxes or form-fields 0, El Ea| ■ Our approach in j^l was different in the regard that ac- 
celeration is possible even if all available eternal fluxes are turned-off: in a sense, the effect of form-fields is replaced 
by a non-trivial "twist" in the internal geometry. Here we considerably extend our knowledge of these solutions by 
presenting new examples, and discuss their properties in a more general setting. 

The rest of the paper is organized as follow. In Sec. II, we review the product space compactification and further 
discuss solutions to the vacuum Einstein equations where the internal product space involves one or more geometric 
twists, or alternatively some non-trivial curved internal spaces. We examine these cosmological solutions broadly into 
three categories: (i) a geometric twist alone; (ii) a non-trivial curvature alone, and (iii) a combination of both. In the 
last case, (iii), we will adopt a different gauge to solve the equations fully. 

In Sec. Ill, we discuss the cosmology on product spaces by using the effective four-dimensional scalars which 
parameterize the radii of the internal spaces and show that when the internal space includes one or more twists along 
with non-flat extra dimensions, then a scalar potential arising from a cosmological compactification possesses a local 
(de Sitter) minimum, whose existence generically predicts a period of accelerated expansion in the four-dimensional 
Einstein conformal frame. In Sec. IV, we discuss the properties of the effective potential in a more general setting, 
which relate our solutions to the solutions given in literature in the context of canonical 4d gravity coupled with 
multiple scalars in an exponential form. In Sec. V, we will present more new examples where the physical 3-space 
expands faster than the internal space. 



The interest in time-dependent supergravity (or S-branes [T8|) solutions arises mainly from the relation of this 
system to string/M-theory and from a 'no-go' theorem E3 which applies if one does not consider cosmological 
solution. In fact, the physical three space dimensions as ordinarily perceived are intrinsically time dependent, and 
if the fundamental description of nature would involve 6 or 7 (hidden) extra dimensions of space, as is required for 
string/M theory to work, then it is reasonable to assume that the internal space is also time dependent. We therefore 
consider model cosmologies described by a general metric ansatz in the following form: 



where the parameters Ti define appropriate curvature radii, and As\ is the metric of the physical large dimensions in 
the form 



II. PRODUCT SPACE COMPACTIFICATION: BASIC EQUATIONS 




(1) 




a(u) 2S du 2 + a{u) 2 ASl 2 k 



(2) 
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so that (f> = <j>(u). Here S is a constant, the choice of which fixes the nature of the time coordinate, u, and dE^. e . 
are the metrics associated with rrij-dimensional Einstein spaces; the values of q = 0, +1, —1 correspond to the flat, 
spherical, or hyperbolic spaces, respectively. 

For the metric Q), the Einstein-Hilbert action is 



/ = 



1G 



-1— f d i+n x^—gR 



1 



16nG 



4+n 



^4 



E 



RlTLi 



(3) 



where ^T. m.j = n. The dots denote, apart from the obvious kinetic terms for the volume moduli 4>i, the other possible 
contributions to the 4d effective potential in the presence of fluxes and/or geometric twists in the extra dimensions. 
When the size of the internal space changes, as in the present case, since 4>i = (f>i(u(t)), the effective Newton constant 
becomes time dependent, in general. This is however not preferable for a model of our universe. To remedy this 
problem, one must choose a Einstein conformal frame in four dimensions by setting |23j 



* = 2 E 



(4) 



The Newton's constant in four dimensions is then time-independent. With Q), the (4+n) dimensional action simplifies 
into the sum of the 4d Einstein-Hilbert action plus an action for the scalar fields (j>i, which determine the size of extra 
dimensions. 



In the product space compactification case, as studied in detail in [23ll24| . we have m t > 2. However, in this paper 
we focus primarily in the case where some of the internal (product) spaces involve one or more geometric twists, so that 
irii = 1 for some i's. For a mathematical simplicity we work in the gauge 5 = 3, in The tt- and ra-components 
of the field equations are 



2K - + 3 — - 3^- 



a 
a" 



- 2k a 4 = 0, 



where prime denotes a derivative with respect to u, k is the spatial curvature and 

m^rrii + 2) A , 2 

i=l 



E 

i>j=l 



(5) 
(6) 

(7) 



In the case of compactification on symmetric spaces taken in direct products, the (scalar) field equations associated 
with the extra dimensions are 



$ + {mi - l) ( j i e- (mi+2)0 '- i: ^ mj ^a 6 = 0, 



where Ui = ti/rf. However, with the introduction of a geometric twist along the internal space, the field equations 
for 6; will be modified. 



To be more precise, let us consider a (4 + n) dimensional metric ansatz, with n = p + 2q + 1, 



d s 4+n — e 



-2$ 



(9) 



where MP is a p-dimensional space of constant curvature of sign t\ = 0, ±1. The remaining 2q + 1 dimensions form 
a twisted product space 



M 



2q+l _ 



= M 2 x M 2 x • •• x M 2 kS 1 , 



q times 



whose metric may be given in the form 

*4 q+ i = r 2 2 e 2 ^ds 2 (M 2 l) + r\e 2 ^ (dz + fw t2 ) 2 , 



(10) 
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with / being the twist parameter. More specifically, 

Q 

ds 2 {M 2 ^ l )= '^(dx 2 + sin 2 Xi dyf), vj +1 = cosxi dy t , 

i=l 

1 I 

ds 2 (M 2 q ) = ^(dx 2 + dy 2 ), w Q = -(x t dy, t - y l dx i ), 

1=1 

8 

ds 2 (M 2q 1 )= ^2(dx 2 + sinh 2 Xi dy 2 ), w_ x = cosha^ dy { . (11) 

i=l 

In the gauge S = 3, the field equations for the scalars, with arbitrary curvature ej, may be given by 

<f>'l + (p - l) (Tl e-( p+2 ^ 1 - 29<#,2 -^a 6 = 0, 

</>2 + oae-^ 1 - 2 ^ 1 ^ 2 -^ 6 - 2j F 2 e^ 1 - 2 ^ +2 ^ 2+ ^a 6 = 0, 

(f)'l + 2< ? F 2 e-^ 1 - 2 (' ?+2 ^ 2+ ^a 6 = 0, (12) 
where a x = ei/rj, a 2 = e 2 /r| and F = (//2)(r 2 /r 2 ). 



A. Compactification on symmetric spaces 



Let us consider the case where cr,; = and F 2 — (i.e. without a geometric twist). The scalar wave equations then 
reduce to fa' — 0. In this case one can easily solve the Eqs. (0,©, with k = 0. First consider a special case where 
fa =const. The corresponding solution is 

cu i , i 2V3 1na fa 

a = a e , fa = const, fa = , — = = ===== = fa H , (13) 

^3z/ 2 + 4q(q + l + " 

where a Q , 0Oi c an< i v are constants. (Of course, the result does not depend on p in this case). 

For 9 > 1, the decompactification of some of the extra dimensions (and compactification of the other) can be a 
slow process. This can be seen also from the solution below, where all harmonic functions <pi (or volume moduli) are 
time-dependent , 

a = a Q e cu , fa — fifa, fa — vfa, fa — ± rjlna, 
up to a (different) shift in fa, where 

(14) 

yJZv 2 + 4q(q + l + v)+p(p + 2)^ 2 + Vi i(y + 2q) 

The constants of integration p, v, rj may be chosen such that each internal space shrinks with time. 

Let us briefly discuss some common features of the above solution by specializing it to a particular model with 
p = 4 and q = 1, so that 

£ 7 = T 4 x T 2 x S 1 . 

As a particular case, let us take fa = — 4(0i + fa) — fa — const. The corresponding solution is then given by 

a(u) = a e cu , fa = ln(a(«)) = -(fa + 0o/4). (15) 

In this example, the 2-space T 2 will grow with time, while the 4-space T 4 will shrink (or vice versa). Unfortunately, 
since the acceleration parameter is always negative, a/a = —2c 2 < (an overdot denotes a derivative with respect 
to cosmic time t, which is defined by t = J a(u) 3 du), there is no acceleration in 4d Einstein conformal-frame. This 
all imply that there is no de Sitter (accelerating) solution in compactifications of a pure supergravity on maximally 
symmetric spaces of zero curvature [17j. However, inclusion of a geometric twist along some of the internal (product) 
spaces, that is crucial to our construction below, circumvents those arguments. 
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B. The solution with a non-zero twist 

Consider the case where <x; = (or a = 0) but / > 0. The equations $F2fy then reduce to 4>'{ = 0, — —q<p2- We 
can then solve the equations (J5J and © simultaneously. The explicit exact solution, with arbitrary values of p and q, 
is given by 

a = a Q (coshx(w - u )) b ° e c o u , 

<t>2 = ~ lncoshx(w — u Q ) + d u = — — , (16) 




where 6 Q = g/4, 



q(3qc + pCj + 4c) _ 3qc + pCj 



d Q = , x= -^^8c d -p{p + 2)c 2 1 -pqcc 1 (17) 

and ao and c are integration constants. Note that not all constants of integration are shown: different (gauge) shifts 
may be taken in (pi ; some of which can be absorbed into the Ti . The value of u is merely a gauge choice and so we 
set u Q = henceforth. The solution l|16f) would be available even if p = 0, in which case as there is no space A4 P , the 
expression for (f>i drops off. 

The 4d cosmic time t is defined by dt — ±a(w) 3 du. The acceleration parameter is given by 

a 5 d = b oX 2 - 2c 2 Q - 46 Q c x tanh X u - & X 2 (1 + 2b Q ) tanh 2 X U, (18) 

where an overdot denotes differentiation w.r.t. cosmic time, t. It follows that solutions will exhibit a period of transient 
acceleration provided that 

2c 2 - 6 (1 + 2b Q ) X 2 > 0. (19) 
Acceleration occurs on the internal w_ < u < u+, where 



-2c Q ±Jx 2 (l + 2b Q )-2cl/b 
tanh X u± = V 77—^-, • (20) 

The number of e-folds during the (transient) period of acceleration, N e = In a(u + ) — In a(uJ) is given by 



N e = -& In 
2 



1 + tanh x u + \ ^ 1 ( 1 — tanh x u x 



(21) 



1 + tanh x u ~ ) \ 1 — tanh x u 

where k = c Q /(xb )- To the parameters given in \17). one has « = ((3q + 4)c + pc 1 )/2%. The number of e-folds iV e 

reaches a maximum when Cj = 0, or cquivalently k = y^ 3 ^ 4 , independently of p. This reflects the fact that the c 1 = 

solution is formally equivalent to the p = solution with no .M p torus and just the twisted (2q+ 1) -dimensional space 
as an internal space. One may show that the maximum number of e-folds is therefore given analytically by 



(^V e ) max = § In i 



3q±4 i / 39+4 

3q 



{g + 2)V3 + 2- y/g(3 g + 4) \ v 3q f (g + 2)V3 + 2 + Vg(3g + 4) 
(9 + 2)^3-2-^(39 + 4) i I (g + 2) V3 - 2 + ^/qj3q~+4) 



(22) 



The number of e-folds increases only marginally when q is increased. 

Let us first briefly discuss some common features of the / > solution by setting q = 2, so that 

b Q = l/2, c = (10c + p Cl )/4 = d + C , X = J 15c " + 3pcc i ~ v ^ + ^ c ?/ 4 - ( 23 ) 
The solution will then exhibit a period of transient acceleration provided that a 2 _c < < ct 2+ c, where 



2p±4y/2p(2p + 5) 
p(3p + 8) 



(24) 
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This implies, e.g. when p — 2, — |c < c± < c. Acceleration occurs on the internal < u < \u + \ where 



, , v — K ± \/2 — K 2 

tanh( X w±) = ^ , (25) 

and k = (10c + pc 1 )/2x- The volume moduli associated with (T 2 x T 2 ) k S 1 are stabilized at late times when 

6c 

VMP + 2) ' 

For p = 2, this yields ci = ±y|c, which lies outside the range that is required for acceleration. In the above example, 

the extra dimensions cannot be stabilized completely, if there is to exist a period of acceleration in the four-dimensional 
Einstein conformal frame. 



C. Solutions with flat and hyperbolic extra dimensions 



In this subsection we focus primarily on solutions that are obtained from compactifications on a product of hy- 
perbolic and flat spaces. In fact, in the F 2 — case, to have a positive potential we need to compactify a higher 
dimensional gravitational theory on spaces that include, at least, one negatively curved factor 

Here we can 

make the problem simpler by considering that (i) the p dimensional space M p is negatively curved while the 2g-space 
M 2q = M 2 x ■ • • x M 2 is Ricci flat (M 2 = T 2 ) (case I) or (ii) MP is Ricci flat and M 2q is negatively curved 
(M 2 = H 2 ) (case II). 

In the first case, the solution is given by 

e c o" 1 
fl(u) %inh X M)V ^^^-^nylnsinhxH, 

4>2 = c 9 u, 6 3 = e, u + 2 In ( ), (26) 

where 6 Q = p/[2(p — 1)] and 

_ p{{p + 2)c + 2qc 2 + c 3 ) _ 3cp + 2qc 2 + c 3 



2 = 1 

X 4p 



4(p-l) ' ^~ 2(p-l) ' 

4gc 2 (3p - 2)(c 3 + c 2 g) + 8q{p - l)c 2 + (7p - 6)c 2 + 3p 2 ((p + 2)c 2 + Aqcc 2 + 2cc 3 )] , (27) 



with c being an arbitrary constant. 

In the second case, the explicit solution is 



e c o u 1 

a ( u ) = Z TT^ - ' 0i= c i u ! 02=c 2 u -^ -InsinhxM, 

(sinhx|M|) z(/ — 1 

fo=c 3 u+2ln( " /2q ~ 1 ) , (28) 



X^ 2 



where b Q = q/(2q — 1) and 



_ q(pc 1 + c 3 + 2(g + l)c) _ pc 1 + c 3 + 6qc 

C °" 2(2q - 1) ' C2= 2(2(7-1) ' 

X 2 = ^ (p [?(3p + 4) - p - 2] c 2 + \2q 2 {q + l)c 2 + \2q 2 {j>c x + c 3 )c + 2(3 9 - l^c.c, + (7q - 3)c 2 ) . (29) 

The solutions given in Ref. [24| may be realized as some limiting cases of the solutions presented above, and fall in 
the category of accelerating cosmologies discussed, e.g. in In the case that W is replaced by S p (or M 2q by 

S 2q ), the function sinh would be replaced by cosh. In this case it is suggestive to study the possible non-perturbative 
(instanton) effects, which may help to decrease the slope of the runway potential. We should however note that, 
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in the case MP — > H p /T, there are no moduli other than the volume modulus, because only modulus of a compact 
hyperbolic Einstein space of dimensions p > 3 is its volume; for a discussion, see, e.g. |28ll34l| and references therein. 



Let us write the scale factor in a canonical form: 



The expansion parameter is 



e c o u 

a(u) = r-. (30) 

(sinh( X M)) h o 



a = ±-^(c - X b coth X u) =±^, (31) 



depending upon the choice of sign in t = ± j a 3 (u)du, while the acceleration parameter is given by 

a = ^ y - b Q x 2 - 2c 2 Q + 4b Q c aX coth X u + b oX 2 {l - 2b Q ) coth 2 X i?J • 



(32) 



The four-dimensional universe is expanding if a > 0. The universe is undergoing accelerated expansion if, in addition, 
a > 0. It follows from l|31fl that a cannot change sign when u < (upper sign) or u > (lower sign), since X b > c Q . 
And the solution will exhibit a period of transient acceleration in 4d Einstein frame provided that 2c 2 + 6 (l — 2b ) X 2 > 
0. Acceleration occurs on the interval u+ < u < U-, where 



2c ± x / X 2(i_26 ) + 2 C 2/6 
coth( xu± ) = . (33) 




FIG. 1: The number of e-folds during acceleration epoch (case I) as a function of the parameter a = c, 2 /c, = c 3 /c, for p = 17, 
q — 2, (d = 26). N e is maximum around (a.,0) — 0. 

The number of e-folds is given analytically by 

a/7 / i 4- 1 \ K+ / rnt.Vi -V7/ — 1 \ K ~ 1 

(34) 



N = ln f f cothxtt+ + / cothxM- - 1 \ K 1 

\ \cothxti- + 1 / \cothxw+ — 1/ J ' 



where 

1 (c, 
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FIG. 2: Tie number of e-folds during acceleration epoch (case II) as a function of the parameter a = c 1 /c, (3 = c 3 /c, for p — 2, 
q = 2, (d = 11). iV e is maximum around (a,/3) = 0. 

Intuitively, the (transient) acceleration occurs far from the cosmological singularity at it = 0. 

To the parameters given in (|27|) or i|29|) . the number of e-folds is of order unity, see Figs. ^ and [21 However, this is 
not a disaster, since the universe would have been expanding prior to the period of transient acceleration, with the 
number of e-folds 



If Am = \u- — u\ < 10~ 36 , then it is not impossible that the total number of e-folds N e + N e > 60, as required to 
explain the flatness problem. In our model, it is possible that the scale factor of our universe remained much larger 
than the size of internal dimensions at the beginning of the accelerating epoch. 



Let us consider the case where, at least, one of the <Ji is non-zero. We also demand that / > 0. In this case we find 
it convenient to choose the gauge 5 = in the metric ansatz J2J), so that u becomes the (proper) cosmic time, i, and 
4>i = 4>i{t). Then upon dimensional reduction from (4 + n) dimensions to four-dimensions, we get 



where G 4 = G 4+ „/Vol(A^"). The kinetic term K is given by 0, after replacing the time derivative ' = d/du by 
' = d/dt, along with the substitutions: mi = p, m 2 = 2? and 7713 = 1. While, the potential term is given by 




(36) 



D. The field equations in S = gauge 




(37) 



y = 7 \ ie -(P+ 2 )0i- 2l ?<fe-03 + A 2 £-P<#>l-2(9+l)<fe-<fe _|_ g-p0i-2(g+2)0 2 +0 3 = ^ _|_ ^ _|_ y p 



(38) 
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where Ai = — e\ p< ^ r ^ , A 2 = — 6 2 -pT and F = (f /2)(r 3 /r 2 ). The analogue Einstein equations are 



1 2 

$ + 3-$-2if-3- = 0, (39) 
a a 

- „ d • /a\ 2 a 2fc . , . 

$ + 3-$ -2 5- =0, 40 

a \ a J a a z 

0'i + 3-0i - -Vi = 0, (41) 
a p 

02 + 3-0 2 - -V 2 - -V F = 0, (42) 

a q q 

03 + 3-03 + 2Vf = 0, (43) 

a 

where $ = p0i/2 + g02 + 03/2. In terms of alternative canonically normalized scalars, which may be defined by 



rJ = \l 2 (q 9 +l) (W>l+2(g+l)02 +03) 



lp(p + 2q + 2) ( 1 



V 2(qr+l) V p + 2q + 2 



the field equations take the following form 



ifi + 3H& + ^- = 0, (45) 



H + K - ka~ 2 = 0, (46) 

along with the Friedmann (constraint) equation 



H 2 = \ {K + V) - ka' 2 , (47) 



where K = \ <^ 2 . The resulting scalar potential is 



y = J \ ie -/3ivi-/32V2 _|_ A 2 e~ 2vi ^' 31 + qF 2 e~ l33Vl+l3iip2 ~ t3::,ip3 , (48) 



where 



A _ ' ' a-v^ItI- < 49) 

9 + 1 P2 VP+29+2 

Clearly the potential has a positive definite minimum with respect to a subset of the ipi directions. By comparing the 
/ = and / > potentials, we see that a geometric twist has a more significant influence on shape of the potential. 
In general, a scalar potential as above, with cross coupling exponents, implies the existence of a local de Sitter region. 
Thus one may expect that accelerating cosmologies with a larger number of e-folds is possible in this case. 

It is rather non-trivial to write the general solution of equations (|45|l - l|47|l with the potential l|48|). in terms of 
proper (cosmic) time t, but any such solutions should be the same as the one from 4 + n dimensional field equations. 
Accepting that the effective potential is a useful tool in this investigation, in later Sections we present various explicit 
solutions that correspond to the above set of equations, but in terms of a new time coordinate, r. 
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E. M-theory phantom cosmology 



Let us first briefly discuss a special case where the scalar potential V(<p) — 0. With Aj = = F 2 , the solution for 
the scalars is 



(50) 



where a dot denotes differentiation with respect to cosmic time t. The factor of is introduced just for a convenience. 
The corresponding kinetic term is 

4 



K = 3C 



a{tr 



where 



C = 6cf + 2c 2 + -C3 + 4cic 2 + c 2 c 3 + 2cic 3 . 



(51) 
(52) 



The effective potential, say V , can have a non-zero contribution if the intrinsic curvature of the physical 3-space is 
non-zero, k ^ 0, see e.g. [23 . I n the k = — 1 case, the solution for the scale factor is implicitly given by 



a(t) 



dx = t + ii. 



(53) 



where £1 is an arbitrary constant. The C = solution is special, as it implies a(t) = t + £ 2 and a = 0. This solution 
critically differentiates between eternally accelerating expansion (if C < 0) and decelerating expansion (C > 0). In 
the case K = V(<p) = 0, the cosmological trajectories can be the null geodesies in a Milne patch (k = —1) of 4d 
Minkowski spacetime, following |3flj 

The choice where all c$ are zero except one corresponds to phantom cosmology, if the non- vanishing Cj is imaginary. 
In this case, since C < (or K < 0), the corresponding solution yields w < — 1, where w = (K — V)/(K + V), and 
also that w— >— last— >oo (for a discussion of phantom cosmology with a single scalar in exponential form, see, e.g. 
Refs. [3(1>II3)- ln our model, however, the choice C < is not physically motivated, because in this case some of the 
extra dimensions behave as time-like, rather than space-like. 



III. MULTIPLE SCALARS AND COSMOLOGY IN FOUR-DIMENSIONS 



In this Section, we give a general discussion on solving scalar field equations with multiple exponential potentials in 
a flat FRW universe. First, we note that for the compactification of classical supergravities on product spaces, with 
or without certain 'twists' in the geometry twist, one can always bring the scalar potential into the form 

s 

V((p) =^A i e"^ E S^ Ai -w- 1< «, (54) 

i=l 

where tpi are canonically normalized 4d scalars. The kinetic term is given by K = IX^V'f- The corresponding 
equations of motion for the scalar fields are 

(pi + 3H<pi + ^- = 0, (55) 



H -TH2--2i (56) 
3 Mi 



while the Friedmann equation is 

p v k 

where a is the scale factor, H = a/a is the Hubble parameter, which represents the universal rate of expansion, 
p v = K + V(ip) is the energy density of the scalar fields, k is the spatial curvature. 

The values of (dilaton) coupling constant Aj are model or compactification scheme dependent, e.g. in the case of 
two scalars arising from hyperbolic compactification of (4+m) dimensional supergravity [23|. we find Ax = 2/A 2 = 

±y/2 + 4/ mi and A 3 = 2^/( mi + m 2 + 2)/[ra 2 (mi + 2)]. Further, A, = - fi m ^-i) ^m_^ ^ SQ that £ . = ^ +1 Qr 

— 1, respectively, for flat, spherical or hyperbolic space. 
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A. The two scalar case 



Consider the potential with 2 scalars (pi,tp2- 

V(<p) = A ie - XlVl + k 2e -^-fi~^f>2 ( 57 ) 

(in units Mp 1 = 1). The exact solution with the single scalar <pi has appeared in j2|| (Ref. provides further 
generalizations, and Ref. j3^,|3!| contain solutions in different time coordinates). In the A2 = case, the field ipi acts 
merely as a non- interacting massless scalar. This case was studied in |38| . Thus we focus here on solutions with more 
than one scalars, where in general Aj 7^ 0. We present solutions with two and three scalars, but the method would be 
equally applicable to higher number of scalar fields. 

In the two scalar case, the late time behavior of the scale factor and the scalars is characterized by 

f (t/t r, fc=o, 

a(t) oc I (58) 
\y/\k\t, ^^0, 



Vl = i (2 lnt - In |^ 0) |) , tp 2 = i Un In t + ^ In \^ \ - ^ In |^ 0) | ) . ( r,!) 1 



where 

_ 2(1 + n 2 ) = A 1 _A 2 

\\ ' M " A 3 
(0) 2(3p-l) /. A 2 \ (0) (3p-l) 2fi 

0\ = — o — 1 - - — U , = : 



7= 72 > M= — \ . (60) 

Ai A3 



^-n^rn 1 -^- ''=ttw (61) 

The late time attractor solutions, with Ai = 0, have appeared before in Ref. 0. Since these asymptotic solutions 
may not be taken to far, below we present a more general class of solutions which are available before the attractors 
are reached. 



B. Accelerating solutions 



With any number of scalars, there would be a period of acceleration (where a > 0, a > 0) before approaching the 
attractor provided that 7 > 1/3. The universe inflates when one of the fields approaches its minimum, where the 
energy density is dominated by the potential energy of scalar fields. In particular, the solution around cpi = and 
ip2 = y~ (c Q is defined below) is given by 



a(t) 



a e 



H t 



0. 



(Vk/Ho) cosh H t, k > 0, 



(62) 



4 (V^k/Ho) sinhffo*, 
where the Hubble rate Hq is characterized by the scale 



k < 0, 



Ho 



7f Y 



(63) 



In a flat universe the solution i|62[) is held only as an intermediate stage in the evolution unless 7 > 1. For k 7^ 0, 
however, the period of acceleration depends on the initial value of tpi, other than the couplings A^. The inflating 
k = — 1 slice that approaches the geodesic a (t) ~ t asymptotically can be eternally accelerating even if 7 < 1, as in 
the single scalar double exponential case [28ll40j] . 

To find the exact solution, with k = 0, one introduces a new logarithmic time variable r, which is defined by 



ln(o(t)), 



-Ai¥>i(t)/2 d £ 
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For k = 0, the field equations (|55(l - (|56|l . with i = 1, 2, reduce to 

1 



^(e + r 1 ) 



pi 



^2 



2^3 
1 



voy/i + /-i 2 (e-r 1 ) 



where a prime denotes differentiation w. r. to r, 

Vb 

and the variable £ satisfies the differential equation 



A 1 + A 2 e c o, 



_,„, A 2 Al-A 2 (Ax-A 2 ) 



111 1 \ 



Ai(Ai - A 2 ) 



f = \A^( 47 ^ 1/2 (U/37 + 1) - (a/37 - !K 2 ) 



(64) 

(65) 
(66) 

(67) 
(68) 



The k = hyperbola 3a' 2 — (1 + /i 2 )(/?' 1 2 /2 = Vb separates the k = +1 and k = — 1 trajectories. Though it is rather 
non-trivial to write the solutions explicitly when k = ±1, one may extract the common physical effects by studying 
the trajectories in a phase portrait, following 0,0]. 

Let us consider the a' > branch, which restricts the solution to expanding cosmologies. For 7 > 1/3, the explicit 
solution is 



;! - +1 tarm^/^T, 



\/37-l 



up to a shift of r about r = 0. Hence 



ao (cosh £t) (sinh (t) 
— In ^(coshCr) 5- (sinner 



1*+ 



4") 



-5-. 



where /3 = \/(l + M 2 )/6, 



1 



3 V37±l' 



C= VVo 



37-1 



4-, 



(69) 

(70) 
(71) 

(72) 



and ao and tpo are integration constants. To obtain the solution for 7 < 1/3, one replaces tanhX by tanX in l|69|l . 
Since the scaling regime of exponential potentials does not depend upon its compactification or mass scales (Aj), ipo 
is actually a free parameter that can, for simplicity, be set to Mp or even to zero. 

The acceleration parameter is given by 



- = 2C 2 e 
a 



[45-6+ - (5- tanhCr) 2 - (6+ cothCr) 2 ] , 



(73) 



which decreases as r increases. This actually implies that the r — > 00 limit corresponds to t — > 00, to the above 
solutions. It follows from l|73|l that, when 7 < 1, the solution will exhibit a period of transient acceleration. This 
occurs on the interval t_ < r < r+, where 



2 _ (2±V3)(V37~1) 



(tanhCr±) _ 

V37+ 1 

For future use, we also note that the shift in cpi during the accelerated epoch is given by 



^r = ~ 



J- In 



cosh (74 



cosh £r_ 

while S(f2 = fJ,Scpi. At late times, r — > 00 (or t — ► 00), we find 



s ln sinhCr + 
+ sinh £t_ 



a 7(7-1) 
— cx 



(3 7 -l) 2 



(74) 



(75) 



(76) 
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An important difference between 7 > 1 and 7=1 solutions is that in the formal case the number of e-folds is arbitrary, 
while in the latter it is fixed, N e = 11.68. Thus, the 7 = 1 solution may be viewed as a transient only since it has got 
a natural entry and exit from inflation; the period of acceleration can be made arbitrary large but not the e-folds! 

The number of e-folds N e during the accelerated epoch is given by 

ftf fT+ I rr+ 

X, = J Hdt = J — dr = J C(^-tanhCT + 5 + cothCr), (77) 



where ti and tj are the starting and ending times of the accelerated expansion. The parameters 5± and £ are constants 
which were defined previously in 1|72[1 . 

For 7 < 1, the number of e-folds during an acceleration epoch is given analytically by 



7 ln (3-\/3)(l + V7) | R 1 ln ^ +1 



3 7 -l ( 3 + V3)(l-V7) V3V37-1 \/3-l 



(78) 



As some specific values, one has N e — 2.62 for 7 = 0.99; N e = 1.5 for 7 = 0.9; N e = 1.16 for 7 = 0.8. It is worth 
noting that the number of e-folds depends only upon the effective value of the acceleration parameter 7, but not on 
the number of scalar fields. 

In the 7 < 1 case, both the lower and the upper limits (t_ and r + ) are fixed in terms of £ and 7, while, in the 
7 > 1 case, only the lower limit of integration (the on-set time of accelerated epoch) is fixed, which is given by 



t_ = — tanh 1 

c 



(2-V3)(V37"l) 



3 7 +l 



(79) 



C. The three scalar case 



It is straightforward to generalize the above results with higher number of scalar fields. Let us consider the potential 
with 3 scalars cpi,ip2,(f3, which is given by 

The solution to the equations 1(55(1 l|56(l . with the potential 1(80(1 . is given by 1(70(1 - 1(71(1 but now 

(81) 



where 



j3 = -= ^/l + ^ + v 2 , V Q = Ai + A 2 e~ c o + A 3 e~ d o, 
V6 



_ A 3 (Ai - A 4 ) - A 5 (Ai - A 2 ) 2 , 2 2 ^ 

v = — , 7= ^ (1 + /T +!/ J 



c = d Q - In 



A3A6 
A3 - v\ 5 



d Q = In 



A2 i/A 3 
A3 A 6 (A 3 - fi\ 2 ) + ^(A 2 A 5 - A3A4) 



AiA 3 



and 



ip 2 = fupi + A 3 l c , (p 3 = v<p\ + A 6 1 d . 

As for the attractor solutions, one reads tpi, ip2 from 159(1 . while ^3 is given by 

2v , 1 1 1 (0)1 , ^5 .(0), A3A4 - A 2 A 5 (o), 
953 = — hit - — ln|^ 3 ; | + — —\n\ip 2 >\ + — — — ln\ip\ '\ 

M *6 A3A6 A1A3A6 



The constants ipf 1 ^ are now given by 



Ai^ 0) 




37-I 


= K 




2 


37-I 


AiA 



A 2 (A^ + A2)- A1A3A5 Uv 



A 3 A| 



A 3 ^ 0) 



A3A6 
2;/ 



A 6 J ' (3 7 - 1) AiA 6 



(82) 
(83) 
(84) 



(85) 
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Notice that only the product A-i(p\ ' is fixed but not each term separately, so the increase of ip\ to ^ can be 

absorbed by rescaling Ai — > Ai (tp^ /tp^). The scalar fields ipi transform under scale transformations in a canonical 
way, i.e. ifi —npi/L. 



IV. DE SITTER VACUUM AND COSMIC ACCELERATION 



A. Exponential potentials with more than one scalar 



Let us first discuss some general features of the cosmological potential (|48|l with F 2 = 0. The corresponding solution 
will then exhibit a period of transient acceleration, with the number of e-folds 



2(n-2) 3 + V3\ATTT-V^T V 3 ^/3(n - 1) - y/3-1 

where n = p + 2q + 1. Thus the number of e-folds increases only marginally when the dimensions of the internal space 
are increased. The late time behavior of the scale factor is characterized by a oc i 7 , where 7 = (n + l)/(n — 1). 

Let us consider some special cases, where one of the scalars ifi takes a (nearly) constant value. In this rather 
restricted case there can exist solutions with a large number of e-folds. One such example is to consider the poten- 
tial (|48|l . with A2 = 0, p = 4, q = 1 and ip\ = b\ =const. Of course this will constrain the evolution of volume moduli 
and hence the 4d effective potential. In this case, the scalar potential is given by 



V = e 



~ bl (Aie -v/ ^ 2 + F 2 e~ 2bl e tp2/ ^' 3V!>/ ^ . (87) 



The explicit exact solution can be found in terms of a new logarithmic time variable t, defined by r 
f dt exp [— (/3 2 (f)/v / 2] ■ The explicit solution is 

a — clq (cosher)" 5- (sinh('r)' 5+ , b, — In— r-r, 



¥>2 = J — In ((coshcY) 5 (sinhCr) 5+ ) , ^3 = </? 2 + — In 



5 V * ' v * ' , , T-» ■ 3 ^9 Af 

where 



5 ± = ^^, C = a/?, (89) 
3(V6±1) V 8 v y 



The total number of e-folds is given by 



.V, = / Ifdf 

-dr+ / —dr. (90) 



T - a' 



The contribution of the first integral will be known once To is chosen, while that of the second term depends on t + . 
As some representative values, with Vq = 1, we find 



N e = 55, 60, 65, • • • , 

respectively, when 8ip2 = 41.91, 45.78, 49.65 • • • (see Fig. (3)). This gives only the minimum number of e-folds, since 
there is an extra contribution from the first integral in (|90|l (i.e. prior to the accelerated epoch), which we have 
dropped here. 

Let us consider some other possibilities, by allowing certain fine tunings among the scalars ipi or the volume moduli 
cj>i. For the potential (|48|l . the late time behavior of the scale factor may be characterized by a(t) oc t 7 , where 
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FIG. 3: The number of e-folds N vs scalar Geld </?2, for the potential \H7\) . 
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g+l , 2(p+2q+2) 
q 9 a (p+2q+3) 
P(?+3)+g+l 

p+2q+3 
_q [_ p+2<?+2 

9+1 



(<p 2 = const, A 2 = 0, Ai ^ 0, F 2 ^ 0) 
(931 = const, A 2 = 0, Ai ^ 0, F 2 ^ 0) 

((^2 = const, Ai = 0, A 2 ^ 0, F 2 ^ 0). 



(91) 



2(<?+l) 2 (p+2q+3) 

For instance, if p = q = 2, the values of the expansion parameter 7, in the above three cases, are 35/18, 13/9 and 
58/81, respectively. In the first two cases, the amount of inflation (or the number of e-folds) can be (arbitrarily) large 
while it is small in the last case. In our model, when all scalar fields (pi vary with time, then a positive potential 
minimum representing the de Sitter phase of our universe can only be metastable. The different choice in the potentials 
could lead to the different asymptotic expansion of our universe. 



B. A model of eternal acceleration 



Let us extend the above discussion by specializing the case with p = 4 and q = 1, but now we take the 4-space Ai 4 
as a direct product between two Ai 2 ,s of arbitrary curvatures. The corresponding lid metric Ansatz is 

ds? x = e- 2 *ds 2 + rJe^da^M* ) + r 2 e 2 * 2 ds 2 (7W 2 2 ) + r 2 e 2ct > 3 {dx 2 + dy 2 ) + r 2 e 2 * 4 (dz + £(xdy - ydx)\ , (92) 

where <£> = <f>i + 2 + 4>3 + <Ai/2, an d ds 2 is the standard FRW metric in the form J2J). Here we work in units where the 
curvature radii n are set to unity, as these variables may often be absorbed in ej (i.e. ej — > Ci/r 2 , where &j = 0,±1) 
and/or in the twist parameter /. In the gauge 6 = 0, upon reductions to four dimensions, the kinetic and potential 
terms are given by 

K = 24> 2 + 20 2 + 2$ + 30 2 /4 + 2^02 + 2<^ 3 + 2^ + fafa + fofa + 3 4 , 

f 2 

y _ _ 6ie -40i-202-203-0 4 _ 4^2— 203-04 _j_ £_ e -20i-202-603+04 ^gg^j 
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In terms of alternative canonically normalized scalars, which may be defined by 

<Pl = 4>1 + 4>2 + 203 + 04/2, f2 = \/3 (&l + ^02 + ^T</>4^ , 

[8/ 1 \ [9 

f3 = \ o 02+O04 , ^4 = -\/o04 ! (94) 



3 

the corresponding field equations reduce to Eqs (J5SJ, (|56|l with ii" = (1/2) ^ <^ 2 . The resulting scalar potential is 

- f 2 

y = _ eie -vi-V3p2 _ e e -vi— v/iva-yfpa _)_ ^g-svi+vT^+vTva— v/fv4 (95) 

4 

For the potential (|95(l . there is a local de Sitter minimum along (f2, which is given by 

V f = - yfa - 2V2^ 3 - f ln ^ A2 ~ 6 W ~ A , (96) 



where 



»2 

f2 p _2 ¥ , 1 -2 v / 6v3 £ = J_ e -6Vl— ^pips— 



ei ei 

Of course, we are allowing here for the possibility that A < 0, B > (i.e. t\ — +1, £2 = —1), in which case there can 
arise two local minima along ip2- With A < (or ^4 > 0) and _B < 0, there is only one local (de Sitter) minimum. 
The solution with A > 0, B > (i.e. e\ — €2 — +1) is, however, unstable. 

Let us consider a special case where ip± = const = 6 . An explicit exact solution can be found in terms of a 
logarithmic time variable r, defined by r = J* dtexp[— ^/|" (^(f)] • The solution is 

a = a Q (cosher) 15- (sinner)' 5- , ifx = b Q , tp% — \/21n [(cosher) 5- (smli^T) - ^] , 

^ 3 = v§ + /5 ln (£f)' ^ 4= vl^ + iM ln ^~ 46 °)' (97) 

where 



S± = —0—, C=l/4^\/l + l/fe-V». (98) 
3(V6±1) V 8 V V 3 v ; 



The late time behavior of the solution is characterized by 

a(t)cct 2 , (pi=bo, if2 = J^hit + const, 

1/33 = In i + const, ^4 = \/21n< + const. (99) 

In terms of the original volume scalars <pi, we find 

e 01 ,e 02 oc f 273 , e 03 cx i" 1/3 , e 04 oc i~ 4/3 . (100) 

In this example the space (T 2 k 5 1 ) shrinks with time, though with different time- varying factors, while the 4-space 
M A expands. 



V. MORE THAN ONE GEOMETRIC TWIST 



A. string/M theory with double twist 



Consider an lid metric ansatz such that 



E 7 = (M 2 ei k S 1 ) x {M 2 2 k S 1 ) x 5 1 . 



(101) 
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We assign to each factor spaces the different time- varying scales: e^ 1 , e^ 2 , • • • , e^ 5 , where, in general, 0i ^ 02, 03 7^ 04- 
We also introduce the two non-trivial twist parameters, j\ and and choose the Einstein conformal frame by setting 
$ = 01 + <fe/2 + h + 4 /2 + fa/2. 

We will begin by presenting a special exact solution where = (so M 2 — > T 2 ) and each 2-space T 2 has the same 
time-varying volume. In this case, in the gauge (5 = 3, the solution will be qualitatively similar to that with a single 
twist, namely, 

a = a cosh(xu) e -^ 8c+c ^ u , 

0x = In cosh \u - -(6c + c 2 )u = 3 , 

h h 
<t> 2 = -0i + c x u + In — + O = 4 + In — , 

Ji h 
05 = c 2 u, fafa = 2 X 2 /a e , (102) 

where 

X = ^48c 2 + 12cc 2 - 8c? - 2c 2 - 4cic 2 . (103) 

This solution will of course exhibit only a period of transient acceleration, with number of e-folds N e = 0(1). However, 
there can arise new solutions with a larger number of e-folds, especially, with a non-zero curvature, e$ =/= 0. 

Let us consider another canonical example, but with six extra dimensions, i.e. Y,q = (T 2 x S* 1 ) x (T 2 X S 1 ). In this 
case, the explicit solution may be given by 



a(u) = a Q cosh(2cu) e 4cu/%/5 , 
0i = 03 = V3cm + lncosh(2cu), 



subject to the constraint 



In |^| + 0o = 04 + hi |^|, (104) 
h h 



8c 2 

/1/2 = — , (105) 



where a Q , c and O are integration constants. This solution can be obtained also from l|102|) after the substitutions: 
c — > — c/v3, ci = and also c 2 = 0, as there is no 05. 



B. Some special solutions 



Let us consider that the internal 7-space is split as in l|101|) . But now we take q + 1 0. Upon reductions to four 
dimensions, we find that the scalar potential due to the non-trivial geometric twists, fa > 0, is given by 



r2 £2 

V(fa,ipi) = h. e -^/\^-\f\v^^^^/\vi-^^o _|_ JJ_ gVlvi— v/fpa-v/fv3— s/Wf 



4 4 
while, the curvature contribution to a scalar potential is given by 



V(ei,VJi) =-eie VsW-vs" V3« V 3^- e2e 
where cp^ are canonically normalized four-dimensional scalars, defined in the following combinations: 



(20! +302 + 203 + 04 + 5 ), <P2 



(20! +203 +404 + 5 ), 



1 / 1 9 

-(01 +303 + 05/2), (^>4 = W -^(801 + 5 ), <^5 = -A/-05- 



(106) 



(107) 



(108) 
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The total scalar potential is then given by V = V(fi, ipi) +V(ei,ipi). 

One may find various exact solutions by specializing the potential V to some particular models, such as 

(pi, ifi2 — const 4>2 = 04 + const, 
f2, <A3 = const => 3 = 4 + const, 
<pi, (f3 — const 302 + 04 — 403 = const. 

Let us consider the first case, where each 2-space Ai 2 will have the same time-varying volume, by further setting 
fi = 0, ei = +1 and €2 = — 1, but f 2 > 0. In terms of a logarithmic time variable r = /* di exp[— \/3 (ps(t)/2] or 
dr = e~^ V3 / 2 6.t, the explicit solution is found to be 

a = a Q (cosher) 5- (sinhCr) 5- , 



ipi = b 1 = const, ip 2 = b 2 = const 

18 
"5 



/18 

</>3 = W — ln[(coshCT) ,5 -(sinhCT)-' 5 +] + </?<,, 



s " = -5|-i71 (1 ° 3 + " l2 - 4h,/2) ' 



where 



(109) 



(110) 



- ^0 f - A / 2i(4-V6) (v ^ fci+fc2)/v ^ 

* 3(^±V3)' ^ V 160 
The late time behavior of the solution may be characterized by 

a(t) cx t 10/9 , y 3 = Int + const, 

V o 

/2 /2~ 

- Int + const, tp 5 = W - In i + const. (HI) 

o v y 

In terms of the original scalars 0,, we find 

e^.e^.e*" oc t 4 / 9 , e 02 , e^ 4 cx t^ 9 . (112) 

In all these solutions the physical three space dimensions expand faster than the remaining ones. In the zero-twist 
case, i.e. /1 = 0, f 2 = 0, we find that, when = —1, the asymptotic behavior of the scale factor is a(t) cx i 2 / 3 . This 
behavior, however, can be different in the presence matter fields, like, dust (u> = 1) and radiation (u> = 4/3). 



C. Slowly rolling moduli 



In this subsection we show that there exist cosmological compactifications of a 4 + n dimensional Einstein gravity 
on twisted spaces of time-dependent metric where 3-space dimensions expand much faster than the remaining n. To 
quantify this, let us consider a decomposition 

£ = (T 2 x T 2 ) x S 1 x (T 2 x T 2 ) x 5 1 , 

by assigning different time- varying scale factors e^ 1 , • • • , e^ 4 , respectively. We also assign the two non-zero twist 
parameters /1 and f 2 , and choose the 4d Einstein conformal- frame by setting $ = 201+02/2 + 203 + 04/2. While it is 
not clear if this model is phenomenologically viable within the string/M-theory context, mainly because the number 
of extra dimensions n is > 7, it nonetheless provides an interesting example where at late times the size of extra 
dimensions can be much smaller compared to the size of the physical universe. 
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Let us begin by presenting a special exact solution where each factor space T 2 has the same volume. The solution 
is most readily written down in the gauge S — 3, which is given by 



a e 2 ( d2 / c+c >eW ' 

1 

'3 ' 



subject to the constraint 



hdu + ln^ + o = 04 + ln^, (113) 
h h 



/ x / a = 2/«o- (H4) 



where ao,c and d are integration constants. This solution exhibit a period of transient acceleration when c 2 > 2d 2 . 
Acceleration occurs on the interval u + < u < w_, where 

1 



u+ = In 

* 24c 



24(2c 2 - d 2 ) ± 24v/3c 2 (c 2 - 2d 2 ) . (115) 



One also notes that 



a = a e( 3c+ -)'" + (-) V4C "e^. (116) 



Clearly, our solutions contain new parameters, the time u = u < at which all space dimensions have comparable 
size, a ~ aoe^. However, at late times (u — > +oo), the ratio between the two scale factors can be such that the size 
of the internal space is much smaller compared to the scale factor of physical 3-space. 

A solution qualitatively similar to the above one arises also for the following decomposition 

E n = (T 2 k S 1 ) x (T 2 x S 1 ) x (T 2 x S* 1 ), 

with the time- varying scale factors e^ 1 , • • • , e^ 6 , respectively. We have been able to find the explicit solution only in 
the case when each 2-space T 2 has the same scale factor, which is given by 



lld 2 \ \ /3 1 



1 12 ' "J eXP U^ e 

1 au (9c 2 -15d 2 )u 

h = 03 = <fc = — e ~ = — , 

a z 2a 

hh 

h = -§\ - du + In — — , 
/1/3 



(XU 



64 = —03 — du + In 



2 



/1/2 

<^6 = ~4>5 - dw + ln— ^— , (117) 

where a = 9c + rf,with c and d being the integration constants. It is possible that the exact solutions that we have 
written above correspond to some special cases, where each 2-space T 2 has the same volume, while, in general, each 
space in the product spaces can have a different time-varying scale factor. In any case, these examples clearly show 
that a universe with 3 + n space dimensions, even of comparable size at early time, could evolve to a universe in which 
the 3 space dimensions become much larger than the remaining n dimensions at late times. 



D. A scalar potential of appropriate slope 



It is well appreciated that a single scalar in exponential form V = Voe~ Xv , with the slope A < v2, can reasonably 
explain the cosmological inflation with the number of e- folds iV e > 11.68 (cf. the discussion in Section III), see, 
e -g-i Whereas such models of string/M theory origin may be quite interesting, no explicit model of this type 

has been constructed so far from compactification (or higher dimensional gravity). Here we give a simple example 
where this may be accomplished. 
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Consider that the internal 7-space £7 is split as 

£7 = M e ei k S 1 

with the scale factors rie* 1 and r2e* 2 , respectively. The parameter w (appeared in JTOJ) may be realized now as / 
(which is a twist parameter) times one-form on j\4 6 . A simplification occurs when j\4 6 is split as M 2 x j\4 2 x M 2 . 
One also chooses the Einstein conformal frame by setting <f> = 30i + 02/2. Let us choose the gauge 6 = 0. The kinetic 
and potential terms are then given by 

1 / -2 3-2 

K = - (240i + -h +60k 

^ = _3^ e -^ 1 -< fc + /!l| e -io^+0 2 _ (n8) 
rj 4 r\ 

In terms of alternative canonically normalized scalars pi, which may be defined by 

(/3l=VT801, ip 2 — ±V6 (01 + 02/2) , 

the kinetic term is K = \{<pi 2 + </?2 2 )- The resulting scalar potential may be given by 

q f2 ^2 

y = __L e -\/2 V i-V 2 /3V2 , e -\/6 V i + V2/3 ^ 



3-. ' ' ( 119 ) 

rf 4 r| 



In the case t\ = +1, we demand / > vl^^Ai) , so that V > when </?i ~ p>2 — 0, but there is no such restriction 
for ei = and t\ = — 1. The potential 1)1191) has a positive definite minimum along <pi (if ei = +1) or along ip? (if 
£i = —1). If any such minimum represents the present de Sitter phase of our universe, then it can be only metastable. 
In this case both the fields (pi and p>2 vary with time. However, a physically different (and perhaps more interesting) 
situation arises when ip2 takes a (nearly) constant value (i.e. 0i ~ —02/2 + const), in which case when the size of S 1 
expands, the size of A4 6 shrinks, providing a "4+1+compact space" type background. 



E. A canonical example 



Finally, as the most canonical example, assume that in the four-dimensional effective theory the scalar potential 
is parameterized by V = Voe~ Xv . One may think of this potential to arise from a compactification of (classical) 
supergravity on spaces with a time-dependent metric p3l |2||, or from a standard flux compactification in string 
theory. The scalar wave equation is then 

i/3 + 3H(p - XV^-^ = 0. (120) 

Let us also include the possibility of a stress tensor for which the mass-energy density evolves as 

p u = -3H(P U + p u ), (121) 

where P u = (lu — l)p u . In general, the energy component p u is not explicitly coupled to the scalar field, but 
gravitationally coupled via the modified Friedmann constraint for the Hubble expansion: 

h2 =U¥ ,+v °^* + ^)- (122) 

Clearly the spatial curvature term k/a 2 acts as a fluid density with uj = 2/3; in this context, a physically relevant 
case is k = — 1 (i.e. cq > 0). The general solution to the above set of equations may be written explicitly in terms of 
cosmic time, t, when lo — 2 (stiff matter), while in other cases, like uj = 1 (dust) or lu = 4/3 (radiation), it might be 
necessary to adopt some other time-coordinates. Fixed point (asymptotic) solutions of the evolution equations with 
uj = 2/3 have recently appeared in Ref. p3|. 

With uj = 2, the explicit exact solution is given by 

1/3 



a I * 6/A 

2 

tp = — In t + In 
A 



t 

X 2 
A 2 V Q 



(> 2 - 6 >/ a2 + a 2 ) 



^^0 + A 2 ^)- 2 , (123) 
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TABLE I: A summary of compactifications, by internal space E n (n — p + 2q + 1), showing the basic characteristics of the 

models. The 2-dimensional space M 2 is either S 2 or H 2 . 



Internal space Number of twists Fine-tuning Compactifying space Decompactifying space Acceleration 
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a linear ref 
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no cosmological event horizon, and so it does not suffer from the problems with the existence of an event horizon 
discussed, e.g., in |44|. The uj — 2 case corresponds to a stiff matter with the equation of state P u = p u . In such a 
medium, the velocity of sound approaches to that of light and hence the cosmology with oj = 2 may differ considerably 
from the usually contemplated scenarios with u> = 1. We will return to the lu 2 case in a future publication. 



VI. DISCUSSION AND CONCLUSION 

In Ref. |17| . Townsend and Wohlfarth studied Kasner type solutions of a pure gravitational theory in 4+n dimensions 
which avoids the no-go theorem of [l6| by introducing negatively curved (hyperbolic) extra dimensions with a time- 
dependent metric. Here we have demonstrated that the same no-go theorem can be circumvented by compactifications 
on Ricci-flat internal spaces. This is possible through the introduction of one or more geometric twists in the internal 
(product) space. That is, a cosmological compactification of classical supergravities on Ricci-flat spaces with a non- 
trivial twist can naturally give rise to an expanding four-dimensional FKW universe undergoing a period of transient 
acceleration in 4d Einstein frame, albeit with e-folding of order unity. 

There was no way in the TW approach to account for a preferential compactification of a certain number of extra 
dimensions. In view of this somewhat discouraging result, one might welcome a modification of this system in which 
some of the extra dimensions shrink with time and/or the runway potential is not too steep. 

We have explained how the scalar potentials arising from compactifications on maximally symmetric spaces are 
modified when there exist one or more non-trivial twists. The combined effects of the non-trivial curved spaces and 
the geometric twists could result in a sum of exponential potentials in the effective 4d theory which involves many 
fields and cross-coupling terms. Such potentials are more likely to have a locally stable minimum with a positive 
vacuum energy. By expressing the field equations in terms of canonically normalized four-dimensional scalars, we 
have obtained a general class of exact cosmological solutions for multi-scalar fields with exponential potentials. We 
also investigated the case where some of the scalar fields take (nearly) constant values (which will then constrain the 
evolution of volume moduli and hence the runway potential in the higher 10 or 11 dimensions) and demonstrated 
that a number of novel features emerged in this scenario, including the possibility of increasing the rate of expansion 
(or the number of e-folds) when there existed a mixture of positive and negative slopes in the potential. The basic 
characteristics of the models (compactifications) studied in this paper are summarized in Table /. 

For all cosmological solutions arising from compactifications on spaces with a time-dependent metric the modulus 
field must give rise to a runaway potential, if it is to allow an expanding four-dimensional FRW universe undergoing 
a period of acceleration. This follows from just accepting the fact that in an effective four-dimensional cosmology 
these scalars can be time-dependent. The rolling of the moduli can be minimized by introducing one or more 
geometric twists along the extra space along with curved internal spaces. So the issue should not be that one needs 
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to find a mechanism to stabilize the moduli completely, but rather that to make any cosmological model arising from 
compactification phenomenologically viable, it is required that the extra dimensions become unobservably small at 
late times by contracting or expanding at a much slower rate than the physically observable dimensions. 

The universe after inflation is dominated by matter fields (radiation and dust). In higher-dimensional supergravity 
theories, supersymmetry requires the presence of specific matter fields, which may hold the clue as to why precisely 
three space dimensions stay large up to the present epoch. In order to fully account for the fate of the extra 
dimensions, it may be necessary to include excitations of the Fermi fields in the cosmological solutions. If so, our 
results are presumably restricted to an era close to the dimensional transition, where the physical three space-like 
dimensions become distinguished from the remaining 6 or 7. 

In this paper we have studied the possibility of generating inflation from higher dimensional gravity on product 
spaces, by introducing only non-trivial curved spaces and/or certain "twists" in the geometry. It would be worthwhile 
to investigate the effects of non-zero (electric) form-fields in spacetime dependent (warped) compactification by turning 
on magnetic field strength (or e quiy alently a geometric twist) and to study their effects on the evolution of the 4d 
spacetime. Some earlier papers |45l |46[ . involving a discussion of the cosmological advantages of string/M theory 
compactifications with twisted internal spaces and fluxes, may be of interest in this context. 

The discussion in 9] shows that the properties of the attractor solutions of exponential potentials (with two scalars) 
may lead to model of quintessence with currently observationally favored equations of state, i.e. w ~ — 1. It may be 
worthwhile to redo this analysis by considering multi-scalar exponential potentials with the slopes predicted by our 
specific examples. One possibility is that desirable slopes can be obtained from the assisted behavior when one or 
more of the fields take a relatively constant value for a sufficiently long time during inflation. An obvious advantage of 
having multiple exponential potentials is that there can exist more than one stage of inflation, some of which happens 
in a vicinity of different minimum of the effective potential. 
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